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Hydrogen  bonds  are  strong  orienting  forces  in  the  organic  solid  state.  A  general  heuristic 
principle  that  has  guided  our  cocrystallization  studies  is  "the  best  hydrogen-bond  donor 
hydrogen  bonds  to  the  best  hydrogen-bond  acceptor".1  Successive  hydrogen-bond 
donors,  barring  severe  conformational  or  crystal  packing  constraints,  hydrogen  bond  to 
available  acceptors  in  rank  order.  Therefore  an  a  priori  knowledge  of  the  relative 
hydrogen-bond  strength  of  the  hydrogen- bonding  groups  within  a  given  cocrystal  should 
allow  for  the  prediction  of  the  hydrogen-bond  connectivity  pattern  in  the  solid  state. 

Within  our  laboratory,  Taft  S  values2  have  been  used  to  assess  the  hydrogen-bond 
basicity  and  pKa  values  have  been  used  to  estimate  the  hydrogen-bond  acidity  of 
compounds  or  specific  funcuonal  groups  within  a  cocrystal  pair.  No  general  correlation 
exists  between  pKa  and  the  free  energy  of  hydrogen-bond  association  across  different 
families  of  hydrogen-bond  donors;  therefore,  the  use  of  pKa  to  assess  the  relative 
hydrogen-bond  acidity  of  donors  is  prone  to  inaccuracy.  Relative  scales  of  hydrogen- 
bond  acidity  such  as  the  ct^  scale3,  which  have  been  derived  from  the  free  energies  of 

hydrogen-bond  association  in  carbon  tetrachloride,  do  exist.  However,  tabulations  of 
are  less  extensive  than  B  and  do  not  include  solid  compounds  of  prime  interest  in 
cocrystallization  studies.  Therefore,  chromatography  was  selected  as  a  potential 
empirical  method  to  assess  the  hydrogen-bond  acidity  of  hydrogen-bond  donor 
compounds  used  in  cocrystallization  studies. 

In  chromatography,  small  differences  in  molecular  interaction  between  a  solute  for  the 
stationary  and  mobile  phases  result  in  macroscopic  differences  in  chromatographic 
retention.  The  logarithm  of  the  chromatographic  capacity  factor  (log  k')  is  proportional 
to  the  free  energy  of  partitioning  of  the  solute  from  the  mobile  phase  into  the  stationary 
phase  and  as  such,  accounts  for  specific  hydrogen-bonding  interactions  and  van  der 
Waal’s  interactions  of  the  solute  with  the  stationary  or  mobile  phase.  In  order  tc  isolate 
the  hydrogen-bonding  contribution  to  the  retention,  two  LC  approaches  were  devised 
that  involve  comparing  the  retention  of  a  given  solute  on  two  matched  stationary  phases 
or  on  two  matched  mobile  phases,  respectively.  The  columns  and  mobile  phases  for  the 
two  different  approaches  are  carefully  selected  to  be  matched  in  their  nonhydrogen- 

1  Etter,  M.C.  J.  Phys.  Chem.  1991, 95, 4601. 

2  Kamlet,  M  J.;  Abboud,  J-L.  M.;  Abraham,  M.H.;  Taft,  R.W.  /.  Org.  Chem.  1983, 48, 2877. 

3  Abraham,  M.H.;  Grellier,  P.L.;  Prior,  D.V.;  Duce,  PP.;  Morris,  J.J.;  Taylor,  P.J.  J.  Chem.  Soc., 

Perkin  Trans.  2  1989, 699. 


bonding  properties  but  different  in  their  hydrogen-bond  accepting  abilities.  As  shown  in 
Scheme  1,  two  polymeric  columns,  a  poly(vinylpyridine-divinylbenzene)  column  and  a 
poly(styrene-divinylbenzene)  column,  along  with  a  90:10  acetonitrile/dimethylsulfoxide 
mobile  phase  are  implemented  for  the  dual  column  approach.  For  the  dual  mobile  phase 
approach,  solutes  are  chromatographed  on  the  poly(vinylpyridine-divinylbenzene) 
column  using  two  binary  mobile  phases,  a  95:5  acetonitrile/triethylphosphate  and  a 
80:20  acetonitrile/triethylphosphate  mobile  phases. 


Scheme  I 
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Both  approaches  were  critically  evaluated  to  see  if  the  difference  in  the  logarithm  of  the 
capacity  factors  (Alog  k'  =  "comparative  retention")  on  the  dual  columns  or  dual  mobile 
phases  was  a  linear  function  solely  of  the  hydrogen-bond  acidity  of  the  solute.  The 
contribution  to  the  measured  Alog  k1  for  each  of  the  various  factors  that  influence 
retention,  such  as  solute  size  (V2/IOO),  dipolarity  (Ttj),  polarizability  ($2),  hydrogen- 

bond  acidity  (0.2)  and  hydrogen-bond  basicity  (62),  was  investigated  by 
chromatographing  selected  series  of  compounds.  The  influence  of  solute  size  on  the 


comparative  retention  was  discerned  by  chromatographing  a  4-alkylaniline  homologous 
series.  The  effect  of  solute  dipolarity  and  solute  hydrogen- bond  basicity  was 
investigated  by  chromatographing  a  series  of  hydrogen-bond  acceptor-only  compounds 
that  covered  a  wide  range  of  dipolarity  and  hydrogen-bond  basicity.  The  effect  of  solute 
polarizability  was  evaluated  by  chromatographing  a  series  of  polycyclic  aromatic 
compounds.  The  dependence  of  the  comparative  retention  on  hydrogen-bond  acidity 
was  tested  by  chromatographing  compounds  ranging  in  acidity  from  aniline  to  p- 
nitrophenol. 

For  each  test  series  of  compounds,  linear  regression  analysis  was  performed  to 
determine  the  multiplicative  coefficient  for  the  parameter  being  probed.  Multiple  linear 
regression  analysis  was  used  to  determine  the  coefficients  for  those  series  of  compounds 
in  which  more  than  one  factor  was  being  probed  and  statistical  tests  on  the  regression 
variance  were  used  to  determine  the  statistical  validity  of  the  effect  of  each  factor  on 
Alog  k\  Values  for  the  various  solute  parameters  -  V2/IOO,  §2.  «2.  and  62  -  were 

taken  from  the  literature.4  Since  the  solute  parameter  sets  each  cover  the  same  range  of 
values,  direct  comparison  of  the  parameter  coefficients  allows  for  the  assessment  of  the 
dependence  of  each  factor  on  the  comparative  retention.  As  shown  in  Table  I,  for  both 
approaches  the  solute  hydrogen-bond  acidity  had  the  greatest  effect  on  the  comparative 
retention.  However,  the  Alog  k'  values  determined  from  the  dual  mobile  phase  approach 
were  less  influenced  by  the  solutes'  nonhydrogen-bonding  attributes  and  therefore, 
provide  a  purer  measure  of  the  solute  hydrogen-bond  acidity.  Principle  component 
analysis  of  the  dual  mobile  phase  Alog  k'  data  for  the  hydrogen-bond  donors  indicates 
that  over  99.9%  of  the  variance  is  explained  by  one  explanatory  variable,  which  from 
the  experimental  results  can  be  assigned  to  the  solute  hydrogen-bond  acidity. 


4  (a)  Bondi,  A.  J.  Phys.  Chem.  1964, 68, 441.  (b)  Li,  J.;  Zhang,  Y.;  Dallas,  AJ.;  Carr,  P.W.  J.  Chrom. 

1991,550, 101-134.  (c)  Kamlet,  M.J.;  Doherty,  R.M.;  Abraham,  M.H.;  Marcus,  V ;  Taft,  R.W.  /, 
Phys.  Chem.  1988,  92, 5244  -  5255.  (d)  Abraham,  M.H.;  Grellier,  Pi.;  Prior,  D.V.;  Duce,  Pi>.; 
Morris,  J.J.;  Taylor,  PJ.  J.  Chem.  Soc.,  Perkin  Trans.  2  1989, 699-711.  (e)  Abraham,  M.H.;  Grellier, 
P.L.;  Prior,  D.V.;  Morris,  JJ.;  Taylor,  P.J.  J.  Chem.  Soc.,  Perkin  Trans.  2  1990, 521-529. 


Table  I 

Dependence  of  the  A  log  k'  Values  on  Various  Solute  Parameters 


Solute 

Matched  Stationary  Phases 

Relative 

Matched  Mobile  Phases 

Relative 

Parameter 

Coefficient 

Value 

Coefficient 

Value 

a  2 

1.6 

1.00 

0.78 

1.00 

f*2 

0.061 

0.039 

t 

t 

* 

*2 

0.19 

0.12 

-0.051 

0.065 

52 

0.61 

0.38 

0.10 

0.13 

V2/100 

-0.78 

0.49 

0.18 

0.23 

f  insignificant  dependence  on  this  variable  as  deemed  by  the  Ehrenson  F  test.5 


Currently,  experiments  are  underway  to  utilize  the  dual  mobile  phase  approach  to  assess 
the  hydrogen-bond  acidity  of  various  hydrogen-bond  donors  used  in  cocrystallizations. 
Specific  attention  will  be  given  to  see  if  the  expected  hydrogen-bond  associations  based 
on  the  solute’s  measured  hydrogen-bond  acidity  occur  in  the  solid  state. 


5  Ehrenson,  S.  J.  Org.  Chem.  1979, 44, 1793- 1797. 
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